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Abstract 

Screening of the topological charge due to he fermion-induced interactions is an impor- 
tant phenomenon, closely related with the resolution of the strong CP and U(l) problems. 
We study the mechanism of such screening in a 'correlated instanton vacuum', as opposed 
to the 'random' one. Both scalar and pseudoscalar gluonic correlators are analyzed by 
means of an observable that minimizes finite size effects. Screening of the topological 
charge is established. This allows us to calculate the r( mass without having to invert the 
Dirac operator. We suggest that this method might be used in lattice QCD calculations 
as well. Our results for the screening of the topological charge are in agreement with the 
chiral Ward identities, and the scalar gluonic correlator satisfies a low energy theorem first 
derived by Novikov et al. |l|. We also propose to evaluate the topological susceptibility 
in the Witten-Veneziano formula not in an infinite box in an world without fermions but 
in an infinitesimal box in a world with fermions. 
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1 Introduction 



Tunneling between topologically distinct sectors of the gauge field, which is described 
semiclassically by instantons p|, ||, is known to be one of the major non-perturbative 
phenomena in QCD. Their importance was significantly clarified in the last few years. In 
particular, a large set of QCD correlation functions was studied [|], || || in the so called 
'random instanton liquid model', RILM. The results are in surprisingly good agreement 
both with experiment (see 0) and lattice data ||. They clearly show that many hadrons 
(including e.g. pions and nucleons) are actually bound by the instanton-induced inter- 
actions, and their properties such as masses (and even wave functions ||) can be well 
reproduced by the simplest instanton-based model. 

Secondly, important new results were obtained from the study of 'cooled' lattice con- 
figurations. ('Cooling' is a procedure that relaxes any given gauge field configuration 
to the closest 'classical component' of the QCD vacuum.) Not only were they found to 
be of the multi- instanton type [|l(|, but recent work of Chu et al. has essentially 



reproduced the key parameters of the 'instanton liquid' picture of the QCD vacuum. In 
particular, it was found in |TT] that the mean instanton separation R « 1.1 fm and the 
typical instanton radius p ~ 0.35/m, about 10% off the values suggested by one of us 



a decade ago |12| . They also demonstrated explicitly how hadronic correlators 'survive 
cooling' and gave quite convincing arguments showing that confinement effects play only 
a relatively minor role. 

In spite of such progress in phenomenology and lattice studies, a consistent quantita- 
tive theory of instanton-induced phenomena is still missing. The model mentioned above, 
RILM, assumes random instanton positions and orientations. It is clearly the simplest 
guess, the first step which ignores all effects related to instanton interactions. It seems 
to be a good approximation for multiple phenomena, but there are also important excep- 
tions. In particular, finite temperature phenomena related to chiral symmetry restoration 
transition cannot be understood without the introduction of strong correlations between 
instantons and anti-instantons |TB[]. Furthermore, in the phase above T c these correlations 
are so strong, that they lead to breaking of the liquid into a set of weakly interacting II 
'molecules'. 

In this paper we discuss some global manifestations of the quark-induced interac- 



tion between instantons in the T = case. It is known that (in the chiral limit) they 
should screen large-scale fluctuations of the topological charge. As it is a very spectacular 
manifestation of dynamical correlations, it is natural to start our new series of studies 
of instanton- induced phenomena by reporting these results first. A detailed report on 
properties of the correlated instanton vacuum, including correlation functions and wave 
functions in particular channels, will be published elsewhere 



Before we proceed, let us remind why fluctuations of the topological charge are of 
interest. It is well known that strong interactions conserve CP to a high accuracy. Two 
ways to understand this 'fine tuning' were suggested. First, the parameter may set 

is 



itself to zero due to the axion mechanism WM. A new suggestions by Schierholtz [16 



that 8 = may happen to be a phase transition point. In any case, there are all kind of 
open questions related to QCD at nonzero values of 6, say whether it is a confining theory 

m 



The second possibility is that the value of 9 may be irrelevant because the total 
topological charge of the vacuum is screened completely. This happens if at least one 
quark flavor is massless, but this does not seem to be the case in the real world. Recently, 



a new proposal was put forward by Samuel flT7| who argued that the interaction between 



instantons may screen the topological charge for non-zero quark masses, below some 
critical mass value. 



Another important aspect is the famous U(l) problem [18, 19], related to properties of 
the 7]' meson. Phenomenologically, its detailed properties are very important, for example 
for understanding of the so called 'proton spin crisis'. Theoretically, a recent discussion 



in connection with the Samuel mechanism can be found e.g. in p0| , |21|| in which the 
exchange of an 7/ was interpreted as a mediator of the Debye-type screening. Generally 
speaking, the 'screening' phenomenon puts significant constraints on the parameters of 
this meson. Finite temperature effects in relation to the interpretation of the rj' particle 
as an inverse screening were discussed in [^2| . 

A famous approach to the U(l) problem is based on the Witten-Veneziano formula 
P3| , P^j relating the parameters of the rj' to the topological susceptibility x- Although it 
clearly explains qualitative features of the phenomenon, we cannot be satisfied by it. It 
is very difficult to make such type of relation quantitative because the two sides of the 



equation 'live in two different worlds' | 20fl . The parameters of the r[ meson are obtained 
from experiment and refer to the real world, while x is calculated either in the large N c 
limit, or in quenched QCD [^5 |, or in other non-screened theories. How one matches the 



unit of energy on both sides is therefore a matter of some convention. 

No convincing lattice measurements of the rf properties in QCD with dynamical 
fermions exist so farQ. In this paper we propose a method that does not require the 
inversion of the Dirac operator and might therefore be useful for lattice simulations. 

The results for the rf correlator in RILM 0, ||, |(| are also unsatisfactory: the cor- 
responding channel has shown to produce too strong qq repulsion. Besides that, the 
random ensemble has obviously no screening, so one cannot get any constraints related 
to it. Therefore, in this paper we study an ensemble of 'interacting instantons', and test 
whether and how the corresponding relations hold. An alternative approach to gluonic 



correlation functions, both in random and correlated instanton vacua, is developed in |28 

We suggest that the best way to clarify the issue is to study the dependence of topo- 
logical susceptibility x(V) defined in the sub-volume V. Its dependence on V, as well as 
on the quark masses mf, are the main objective of our work. The screening of the topo- 
logical charge implies that, for large V, x(V) ls proportional to the area of V instead of 
xiy) ~ y i (see also [[26[ where this phenomenon was discussed for the Schwinger model). 



However, for small volumes one always has that x(Y) ~ V • The cross-over point between 
these two regimes defines a screening length, identified as the inverse rf mass. 

In section 2 we define the instanton liquid model and discuss the parameters used 
in our calculations. In section 3 the Debije cloud of a topological charge is studied and 
in section 4 the volume dependence of the topological susceptibility is evaluated. These 
results are analyzed in the context of correlators that follow from effective field theory 
(see appendix A). In section 5 we discuss the scalar gluonic correlation function which is 
related to the fluctuations of the total number of instantons in a given volume. Concluding 
remarks are made in section 6. 



lr The best one so far is [B7L where the r( mass was calculated in quenched lattice gauge theory using 
a rather indirect method. 



4 



2 The model 



The basis of this work is the 'interacting instanton approximation' (HA), formulated as a 
particular statistical model amenable for numerical simulations |29], |30], |3]J or analytical 



studies [32], ^]). The partition function for N/2 instantons and N/2 anti-instantons is 
generally given by 

N N f 

Z = [f[ dtonlzidpinip!) it det(T + im f ) exp[-/3(pA) £ Sff], (2.1) 

J 1=1 f KJ 

where fti,zi and pj, denote the orientation, position and the size of pseudoparticle /, 
respectively. The size distribution p(p) contains a factor ~ p 6 ~ 5 (with the standard 
beta-function parameter b = yAf c — fiVy) due to Jacobian of the transformation to 
collective coordinates and the leading quantum corrections (/3(pA) = —b\og(pA). The 
QCD parameter A will be determined phenomenologically. 

The main dynamical effects we are going to study are due to the determinant of the 
Dirac operator. In the subspace of the zero modes it can be expressed in terms of the 
overlap matrix elements T, which for the 'streamline' gauge fields [33J] are given by [EJH] 



T AI (p lPA ^ = {p lPA f* J d 4 x^l A {x)iD^ {x) = l -Tr (t+R^U^Ua) F(X), (2.2) 
where the scalar function F(X) is defined by 



oo 



F(A)=6 io ( r + l/A)3/2( r + A)5/2- ( 2 - 3 ) 
It depends on the conformally invariant combination A 

A = ^ 2 + - + -VU(^ + - + -) 2 - 4 ^ (2-4) 
2 V Pa Pi J 2 \ Pa pi J 

where R = Rj ^/pYpX- Asymptotically, for large R/p, the Tai ~ 4p 3 /i? 3 (the geometric 
average of pi and pa id denoted by p). 

Considering the interaction due to gauge fields, let us first mention that we deviate 
from previous works^ by using the interaction derived in [|34| from the so called 'streamline 
equation' |36| . They correspond to the true bottom of the valley for the I /-configurations. 



2 They used variational trial functions for the gauge fields known as the 'sum' or 'ratio' ansatz. 



This interaction differs from the previously used ones in one important aspect: for one 
relative orientation the repulsive 'core' at small distances is gone. Our numerical simula- 
tions have confirmed earlier suspicions that it leads to an 'overcorrelated liquid' of close 
II pairs. It is still possible, that one may get a phenomenologically acceptable ensemble 
by a significant increase in the density of instantons. However, to spend a major por- 
tion of computer time simulating pairs that are not even semiclassical fluctuations is not 
practically possible. 

Thus, it seems that the original hopes to stabilize the ensemble at the purely classical 
level []32| are not fulfilled. Presumably quantum effects (especially subtraction of per- 
turbative contributions, relevant for close instanton-anti-instanton pairs with a strongly 
attractive interaction) will generate the effective repulsion, and two recent developments 
should be mentioned in this context. First, Diakonov and Petrov calculated the instanton 
interaction from the semiclassical quark scattering amplitude [37] and found a short range 
repulsion. However, at the moment it is not yet clear to what extent this interaction is 



of relevance to the present problem. Second, some quantum corrections, discussed in [p8 
using the scale anomaly relation, also lead to some effective repulsion. 

Lacking a detailed understanding of these effects, we have introduced a phenomeno- 
logical repulsive cor^\ to the streamline action S SL (X) 



A 

A 



S% = S sl (\) + —Tt1 2 U F 1 U a 1 2 U a 1 U i , (2.5) 



STr = ^UUT'U'^U^Uj, (2.6) 

for pseudoparticles of the opposite and identical kind, respectively. (In the latter case the 
classical interaction is absent and only the hard core is included.) The strength of the 
core is denoted by a free parameter A, and the color traces are needed to to reduce the 
magnitude of the interaction if two pseudoparticles belong to different SU(2) subgroups 
of SU(N C ). 

In principle, A determines the total instanton density in terms of A.q CD . Unfortunately, 
it is not known with sufficient accuracy, and therefore it is logical to do the same thing 
as on the lattice, namely tune those two parameters to some measured masses, and then 



3 Let us point out an analogy to the famous problem of nuclear matter saturation: although there 
is no doubt for the existence of a repulsive core for NN forces and its decisive role in the problem, its 
properties and physical nature is debated even today, after decades of investigations! 
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consider others as predictions. This strategy will be used in a subsequent work on hadron 
spectroscopy in the IIA. 

However, in this paper, related to only a limited number of issues, we adopt a simpler 
normalization prescription. For all masses, number of colors or flavors we (i) fix the total 
density of the instantons at n — Ifm 4 , roughly corresponding to the standard gluon 
condensate value and to lattice measurements; (ii) the repulsive core strength is fixed to 
be A = 128, for which we obtain a quark condensate that is close to the empirical value. 
(Note that in this region the dependence on A is relatively weak.) 

We simulate the partition function with a Metropolis algorithm for an ensemble of 
N/2 = 32 instantons and and an equal number of anti-instantons in a box of size 
(2.38/m) 3 x 4.76/ra. We average over typically 1,000 statistically independent config- 
urations. 

3 Screening of the topological charge 

Let us start with the evaluation of the topological charge 



in a subvolume ('slice' of the box) H(l 4 ) = L 3 x / 4 . It is not an easy task on a lattice, but 
for our instanton ensemble this simplifies to just counting instantons and anti-instantons 
with centers in the subvolume. 

Our aim is to study the details of the screening phenomenon such as the size and shape 
of the corresponding 'Debye cloud'. Let us take the center of one of the anti-instantons 
as the origin of the coordinate system and measure the charge distribution in the rest of 
the system. In Fig. 1 we show results of such studies for N c = 3 and Nf — 2 (both u,d 
quarks with equal masses, given in the label of the figures). Since the total charge in the 
box is set to zero, the 'compensating charge' contained in the slice H(l 4 ) for large enough 
L 4 should approach 1 (the anti-instanton in the center is not counted). 

Our data (points) do show a strong screening phenomenon, namely that most the 
compensating charge is contained in rather small slice of the box, with a width of only 




(3.1) 
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about 1/3 fm. That means that the correlations between instantons and anti-instantons, 
responsible for screening, are actually very short range ones. The relevant hadronic ex- 
citation, the pseudo-scalar flavor singlet channel, should therefore be relatively heavy. 
Thus, our ensemble solves the U(l) problem at least qualitatively. 

One more qualitative phenomenon to mention here is the observed charge oscillations 
around the central charge. Those are well seen for m = 0, and fade away for larger 
masses. The oscillations are clearly non-statistical, and their appearance suggests that 
in the 'instanton liquid' there exist correlated clusters of various sizes. Let us recall in 
this connection, that in many ordinary liquids one finds local order characteristic to a 
crystalline phase, and that clusters can be as big as containing hundreds of atoms. This 
interesting topic clearly deserves further studies. 

Doing a more quantitative analysis, we have fitted these data by the function 

Q(Z 4 ) = a 2 (l - exp(-m/ 4 )) + (1 - a 2 ) y (3.2) 

The fraction a is the magnitude of the 'compensating charge', if a = 1 the charge is 
screened completely. The second term corresponds to the remaining charge which is 
assumed to be distributed homogeneously over the box. The first term is obtained by 
integrating the correlation function (|A.13| ) for fixed — = I4. In the effective low 
energy field theory (discussed in appendix A) the mass m is identified as the rf mass. 

The obtained values for a are 1.0, 0.99, 0.96 and 0.80 for the 4 masses used (m = 
0.0, 0.1, 0.2 and 0.4 in units of fm' 1 ), with statistical errors about 0.05. Thus, the 
non-screened charge (1 — a 2 ) ~ m 5 with the power 5 = 1 to 2. This conclusion should 
be compared to the conventional wisdom, according to which it should be proportional 
to the lightest quark mass. We do not know whether there is a discrepancy, or that it is 
due to finite size effects (the 1/L term in ( |3.2|) ). In the scenario proposed by Samuel, the 
effect should vanish at a non-zero critical mass (1 — a 2 ) ~ (m — m c ) 5 . Our results lead to 
an upper limit of about m c < 20MeV. 

The value of the second parameter, equals 673, 623, 555 and 518 MeV, in these 
4 cases, respectively^. The result is that singlet mass decreases with the current quark 

4 We remind the reader that the 'MeV makes sense in those un-physical theories only after we have 
defined them, rather arbitrarily, by keeping the density of instantons fixed at some value. Comparison 
to lattice simulations is possible, but one has to keep the same convention. Th physical situation will be 
discussed below. 
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mass and suggests that the screening length becomes larger if one is increasing all quark 
masses. 



4 Fluctuations of the topological charge in a subvol- 
ume 



In this section we switch to more conventional observables, which have been discussed in 
both theoretical and lattice literature. They are closely related to what was done above, 
but at the same time they do not demand the charge to be localized. 

Although the average value of the topological charge in some volume is, of course, 
zero, one may study its fluctuations {Q 2 )- More specifically, we study the dependence on 
the the finite volume V4, its surface A3, etc., 

(Q 2 ) v = C v (m) V, + C A (m) A 3 . (4.1) 

(The subscripts refer to the dimensions of the manifolds). 

We have already mentioned, that the attention in literature has been focused on the 
large volume limit of that quantity, the topological susceptibility 

X = lim (4.2) 

V4-»00 1/4 

The dependence of x 011 the quark masses m is the same as for unscreened charge 1 — a 2 
discussed above: the standard scenario with the 'screening of the topological charge' 
implies that x ~ m at small m, while in Samuel's scenario x = for all quark masses 
below some critical value m < m c . 

Our idea is to study not only the limiting value of the topological susceptibility but its 
dependence on the subvolume. The particular behavior to be observed is very informative, 
and can be used to extract additional information about vacuum properties (and hadronic 
spectroscopy). Specifically, as we try to separate volume from surface effects, we take a 
segment of the hypertorus of variable length Z 4 along the 4-axis, i.e. the slice H(l A ) 
discussed in previous section. Its total volume L 3 / 4 is proportional to Z 4 , while its surface 
area L 3 is independent of it. Because our box size in the 4-direction is L 4 = 2L, the 
effective range of / 4 is < / 4 < L. 



Furthermore, as we have seen in appendix A, the topological correlation function 



Hp(x-y) (see (|A.13|) ) can be connected with parameters of the SU(3)-scalar pseudoscalar 
mesons, the famous ?/— particle and the corresponding excited states. 

The fluctuations of the charge Q(Ia) in the box H(U) are obtained by integrating the 
correlation function (|A.13| ): 

ph 
-h 



K P (IA = (Q(h) 2 ) = L 3 [ d 3 x t dtih - \t\)Iip((x 2 + 1 2 )-* 

JLfi J-U 



(4.3) 



This integral vanishes for Z 4 — > oo if any quark becomes masslessQ (see eq. (A. 15) in 
appendix A). For finite U this integral can be done as well; for sufficiently large L ^> 
1/m the range of the spatial integration can be extended to infinity with high accuracy 
0(exp(— Lm^')) and the result is 



(Q(h) 2 ) = L 3 n 



U + cos 2 



m 2 fl — exp^—m^U) 



L + sin d> — - 



mf ( 1 — exp(— rrirjU) 



-U 



rat \ m r 



For I4 — > we find 

(Q(h) 2 } = nL%. (4.5) 

One can easily convince oneself that this result is independent of the specific shape of the 
sub volume, and in general we have 

h M = » (4.6) 



Together with eq. (|A.14|) , this results in the following formula for the 7/ mass 



2^2 2 2iV, {Q{Vf) 
m v , + m v , - m K = hm , (4.7) 

which is our analog of Witten-Veneziano formula, with the topological susceptibility re- 
placed by local fluctuations of the topological charge. It should hold for volumes in 
the window between that of a single instanton 0(p 4 ) and the inverse of their density 
0(l/n) = OCR 4 ). 

In terms of the flavor singlet spectral function, the correlation function ( [4.4p implies 
that neither the contribution of the 1 nor the contribution of the continuum have been 



3 A similar cancellation takes place in many cases, e.g. in (2-d) Schwinger model (see |26| for details). 
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taken into account. We only calculate the non-perturbative contribution to the correlation 
function Up P = Hp — H p ert ' which implies that the contribution of the continuum and the 
resonances appear with opposite signs in Up P . In particular we will find a cancellation 
between the i and the continuum. In the fits to be discussed below we included the 
contribution of the continuum, and in all cases we found that the best fit was obtained 
for zero coupling of the continuum. From now on we do not include it in our discussion. 

Our results for the fluctuations of the topological charge, Kp(U), in a slice with length 
Z4 are shown in Figs. 2, 3 and 4. The first qualitative point we would like to make is 
that the data points clearly show the expected transition from a linear dependence on Z 4 
(for small box size) to a constant for larger boxes (in case the quark mass is zero). These 
data point deviate drastically from the full line corresponding to randomly positioned 
instantonsj], the behaviour that holds for the RILM used in 0, [|, || . Thus, the topological 
charge is in fact screened and its fluctuations become a surface rather then a volume 
phenomenon, as soon as the quark-induced interactions are included in the statistical 
sum. 

In the upper figure of Fig. 2 we compare the cases with one and two massless quarks, 
for N c = 2. For massless quarks the correlation function (4.4) reduces to 



(OWI^ = 1 - eXP( -" V ' 4> . (4.* 



The dashed curves show a fit of this function to the data points. For Nf = 1 and Nf = 2 
we find m'^ = 535MeV and = 756MeV, respectively. The ratio of the masses is 1.41 
which is right on top of the theoretical expectation that m' ~ yNf (see ( A.14j) ). If we 



study the susceptibility (Q(h) 2 ) as a function of the quark mass, which effectively reduces 
the number of flavors, we therefore expect on the one hand that the rf mass decreases, 
whereas, on the other hand, it increases because of the explicit contribution of the quark 
masses. 

The quark mass dependence of the topological correlator is studied in Fig. 3. In this 
case the number of colors is 3 and we consider two flavors with equal mass (same instanton 
configurations as in Fig. 1). For the quark mass we refer to the label of the figure. 



6 The functional form of this curve, which follows immediately from the binomial distribution of the 
instantons, is given by NQ^/L^)^ — (I4/L4)). 
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In the case of equal quark masses, the fitting formula (that follows from eq. |A.13|) is 

(Q(kf) = "i 3 U(l " + (1 " a 2 ) 1 - CXP( -" V ' 4) ) - (4-9) 

The value of a obtained from a least square fit is 0.0, 0.0, 0.0, 0.19 and equals 580, 
520, 430 and 340 MeV for quark masses of 0.0, 20.0, 40.0 and 80.0 MeV, in this order. 
As we already observed in previous section, our results set a limit on the existence of a 
critical mass below which we have complete screening. Again with the disclaimer that the 
individual parameters are not well determined by the fit. Only the height of the curves is 
well-determined. Our values for the r/ mass are somewhat less than in previous section, 
but show the same trend that it decreases for not too large quark masses. 

Finally, in Fig. 4, we show results for realistic quark massesf], m u = = lOMeV 
and m s = 150MeV. The dashed line corresponds to effective theory with realistic meson 
masses and t] — ^'-mixing angle. The full line is a fit of eq. ( [4.4|) to the data. Both the pion 
and the kaon mass are found to be zero, whereas = 605MeV. Also in this case the 
individual parameters are not well determined by the fit. In particular, the fitted curve 
is very insensitive to the value of the strange quark mass. In fact, the best fit misses the 
strange quark contribution to the . Note that in naive chiral perturbation theory this 
meson should possess a strangeness-related mass m s ^ ange = m. n /y/2 = 388 MeV. The sum 
of these two numbers are right in the correct region. 



5 Compressibility of the instanton liquid 



Let us now proceed to studies of a related gluonic correlation function 

2 2 

Ks(x - y) =< JL- 2 FF(x)^L- 2 FF(y) >, (5.1) 

containing the action rather than the topological charge, and therefore counting instantons 
and anti-instantons equally. The total topological charge in a subvolume V is given by 

2 

N(V) = [ d 4 x^FF(x), (5.2) 

7 Of course, for any finite volume simulation the value of to/ cannot be taken arbitrary small: the 
minimum possible value of m is the point where the to— derivative of the condensate vanishes. For 
feasible volumes we work with mj = 0.05 « lOMeV, being small compared to hadronic masses but still 
considerably larger than the light quark masses in the real world. 
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which satisfies that (N(V)) = nV . For pointlike pseudoparticles we have 

.2 N 



g 



FF{x)=Y J 5 A {x-X I ). (5.3) 



32tt , , 



For small \x — y\ only terms involving the same instanton contribute and we find thatQ 



N-six — y) = n5 4 (x — y) for \x — y\ — > 0. (5.4) 

In reality the instantons are not pointlike, and show a repulsive interaction. This leads 
to the low-energy theorem |], [32], |4[] 

([N(V)} 2 } - (N(V)} 2 = ±(N(V)), (5.5) 

for V —>■ oo. The compressibility of the instanton liquid is thus equal to 4/6. As can be 
seen from (|5.4|) , for small volumes the factor 4/6 is absent. How can we understand this 
in terms of the physical states contributing to the correlation function? The lowest lying 
intermediate state is a 'scalar glueball', which we will generically call o (not to be mixed 
with the a of sigma-model, at 600 MeV). Using a somewhat different convention than in 
the pseudoscalar channel the relevant matrix element is 

(0\jL F F\a) = \ a . (5.6) 

Including the contact term and the contribution of the continuum with coupling constant 
A c , this results in the correlator 

POO 

n s (x-y) = n5\x - y) + \lD(m a) \x - y\) + \ 2 C / dsD(^, \x - y\)s 2 . (5.7) 

J SO 

This result contains both perturbative and nonperturbative contributions. However, 
the correlation function that is obtained from the instanton liquid calculation only receives 
contributions from the latter sector. Therefore, we should compare out results to n^ p = 

— n^ ert '. The spectral representation of the perturbative contribution is given by 

POO 

m crt - = X 2 / dsD{^, \x - y\)s 2 . (5.8) 
Jo 



Note that in [28] only the correlation function at finite distances x is considered, so the local term 



we discuss now is not seen. However, we need it because we consider the integrated correlator. 
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In the same way as for the topological correlator we define the correlator, Ks{h), as 
the variance of the total field strength inside the box H(l^). Numerically, this correlator 
can be evaluated again by simply counting the number of instantons inside the box. It is 
related to II<j P by 

rh 
-h 

Using the the correlators ( |5.7| , |5.8|) we find 



K s (h) = L 3 f d 3 x ( h dt(l A - |t|)Il£ p ((£ 2 + t 2 )5). (5.9) 



K s (h) = nL% + \lL\ k - 1 ~ ex P(-^) ) _ xlL s r (k _ LggH^ )^, 

m a Jo \/s 



(5.10) 

For large volumes the low-energy theorem fl5.5| ) should be obeyed. This leads to the 
relation 

n + \l-\ 2 c = ^n, (5.11) 

which shows that in the scalar case (in contrast to pseudoscalar one treated above) it is 
essential to include the contribution of the continuum. (Recall also that in the pseudo- 
scalar case the sign of contribution from the rf resonance is negative due to the anti- 
hermiticity of FF.) 

In the lower figure of Figs. 2 and 4 and Fig. 5, we show results for the scalar correlators. 
The full lines in all figures correspond to a random positioning of the instantons. For 
infinite volume we would get a straight line, but finite size corrections turn it into a 
parabola. The numerical results are represented by dots, and the dashed line is a fit using 
the function 



+ -|-(l + A 2 -^)F(/ 4V ^ )^, (5.12) 



2L \b 2L J m a 

where the 'missing continuum function' F follows from the integral in (5.10) 

F(x) = -\ + ^~^ + i + ^) exp(-x)). (5.13) 

This fitting function both reproduces the short distance behavior ( |5.4j ) and the low-energy 
theorem (|5.5|). 

Results for N c = 2 and Nf = 1 or Nf = 2 are shown in Fig. 4. Two fits with 
comparable x 2 but very different parameters are found. First, we a fit with A = and 
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y/s = 1110 MeV for Nf = 1 and 715 MeV for Nf = 2, respectively. In this case mo- 
is indeterminate. A second minimum is found for A = 1.22 and m a = y/s Q = 2.4 GeV 
for Nf = 1 and A = 1.32 and m a = v^o = 1-9 GeV. In both cases, the bulk of the U 
behavior is determined by the low-energy theorem, and the resonance and tcontinuum 
only contribute significantly for small volumes. 

In Fig. 5 we study the current quark mass dependence for N c = 3 and two flavors of 
equal mass. Again we find two different fits. First, for current quark masses equal to 0, 
20, 40 and 80 MeV we find a fit with A = 0.0 and y/s Q equal to 1160, 1295, 1245 and 1525 
MeV, respectively, and, second, we find A = 1.56, 1.60, 1.44 and 1.72 and m a = ^/so = 
2.6, 2.9, 2.6 and 3.7 GeV, in the same order. In both cases, the bulk of the correlation 
function is determined by the low-energy theorem ( |5.5|) . 

Finally, in the lower figure of Fig. 3, we show Kg for realistic values of the quark 
masses (see caption). As a best fit with A = we find that y/s Q = 1320 MeV. For the fit 
with m a = i/Io we find a mass of 2.97 GeV and a coupling constant of A = 1.53. These 
results are consistent with the factor 4/6 in the low-energy theorem. 

We want to emphasize that the value A = and the equality m a = ^/sq are a result of 
the fit. The fits quoted here are not the only fits; there are other minima of x 2 which we 
dismissed as unphysical. For example, in the last case is also fitted by ^/so = |m CT = 200 
MeV and A = 4.5. For this ratio of the masses the asymptotics of the fitting function 
becomes independent of A. Therefore the conclusion that the bulk of the results are 
described by the low-energy theorem remains unaffected. 

6 Conclusions 

We have presented a detailed study of the volume dependence of the fluctuations of the 
topological charge in an Interaction Instanton Approximation to the QCD vacuum. As 
was expected, we found that these fluctuations are screened for a sufficiently large box and 
sufficiently small quark masses due to the fermion-induced interaction between instantons. 
The phenomenon disappeared gradually as the quark mass increases from zero. An upper 
limit for the Samuel's critical mass value is set to about 20 MeV. Furthermore, the results 
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for a correlated instanton liquid have shown that correlations are essential for reproducing 
the m' v and its dependence on the number of colors and flavors in the theory. 

As a by-product, we have obtained a formula similar to the Witten-Veneziano formula, 
but the topological susceptibility that enters into it is not the large volume limit but rather 
the small volume limit in the presence of dynamical fermions. 

We also present results for the fluctuations of total number of instantons in the finite 
box. The results clearly indicate that the instanton liquid is not an ideal gas, but a liquid 
with an important repulsive interaction. The results agree very well with the NSVZ low- 
energy theorem, and provide constraints for the parameters of the physical particles in 
the scalar channel. 

Finally, let us emphasize that the method to evaluate the r( parameters proposed in 
this work can be applied to lattice QCD calculations with dynamical fermions as well. 
Straightforward evaluation of the quark propagators can be substituted by measurements 
of the topological charge in a finite box, which may be much easier to do. 
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Appendix A: Topological correlators in effective field theory 
The topological charge density is defined by 

Q(x) = -^FF{x), (A.l) 



where F is the field strength tensor. The topological charge Q(V) in 4-volume V is 
obtained by integration of Q(x) over V. The fluctuations of the topological charge in a 
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given box are determined by the topological correlation function defined by 

ri FF(x)). 



U P {x-y) = (^FF(x) !, ~ 



(A.2) 



32vr 2 " "~ / 32n 2 

In this appendix we will derive this correlator from the effective low-energy chiral La- 
grangian [RjJ and its interaction with the topological charge Q. This derivation is based 



on work in references [19, |20], 5T|. The main additional ingredient is that we include the 
77 — 77' mixing in our derivation. 

The topological charge density couples to the pseudoscalar singlet meson channel, and, 
for a diagonal mass matrix with m u = m^, only the 77 — 7/ mixing part of the effective 
pseudoscalar Lagrangian is of relevance^. This part, which involves only the diagonal 
SU(3)— flavor fields O an d 08, is given by [|4], |3^, 40], |33], 41 



off 



w 2N f 



0, 



with 



£(00, 08 



1 

2 



and the square of the mass matrix |24 
M 2 = 



-rr? K + \m\ 



m. 



V 2 + M- 



4 2 

1m 1 



m 



2 

3777,2 



(A.3) 



(A.4) 



(A.5) 



w ■$ k 

The normalization of the fields and the coupling constant / will be discussed below. 
In the dilute gas approximation the instanton partition function is given by f42l , |19 



E 



K N++N- N 



NJN_ 



PJ d 4 ^ exp(- / d 4 xC eS ). 



(A.6) 



The constant k was calculated to one loop order in ||. Treating the instantons as pointlike 
objects, the topological charge density is given by Q(x) = J2Qi<5(% — This allows 
us to perform the sum over N + and N_ in the partition function. Our final effective 
Lagrangian is obtained by expanding the resulting cosine-function to second order (weak 
field approximation): 



C 



2N 

, [f = -2Vk + k-^/0o + £(0 O , ^ 



(A.7) 



9 In fact, there is another state with the quantum numbers of the rj' , 7y(1430) and formerly called t, 
which gives a comparable contribution to this correlation function (see discussion in Qj). We do not 
include it here, because it is hardly possible to extract any information about this state from the data 
sample discussed in this paper. 
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The constant k can be traded for the average pseudoparticle density by differentiating 



both (|A.6|) and ( |A.7|) with respect to k. We find 



2k = (y) ee n. (A.8) 

The topological correlator (|A.2|) can be obtained by integrating out the meson fields. 
However, it is a somewhat simpler to write the correlator in terms of functional deriva- 
tives with respect to 0o of the first term in the Lagrangian ( |A.3| ) and performing the 
differentiation after the sum over the instantons, i.e. on the second term in ( |A.7| ). The 
correlator can be expressed as 

2 N 

(Q(x)Q(y)) = n5\x - y) - n 2 -^(0 o (a;)0 o (y)>. (A.9) 

The latter expectation value is most conveniently obtained by diagonalizing the mass 
matrix A4 (including the topological contribution) with eigenvalues identified as m n > and 
m v and eigenvectors expressed in the r/ — r\ mixing angle according to 

\rf > = cos 0|0o >+ sin 0|08 >, (A.10) 
\r] > = cos 0|08 >— sin 0|0o > . (A. 11) 

For the singlet correlation function we obtain 

(M x )Mv)) = cos 2 ((f)) D(m v ,\x - y\) + sin 2 (<j))D(m v \x - y\) (A.12) 

resulting in the topological correlator 



(Q(x)Q(y)) = 7^ m top $ ( x -y)~ ™% op (cos 2 (0) D{m Tf , \x - y\) + sin 2 (0)D(m^, \x - y\)) 



(A 



where the scalar propagator is given by D(mx) = mKi(mx) / (Air 2 x). For convenience we 
have introduced the topological mass 



m t 2 o P = = m 2 , + m 2 - 277-4, (A. 14) 



where the latter identity follows from the invariance of the trace of the mass matrix 
under diagonalization. This correlation function can also be derived from the pseudoscalar 



singlet spectral function with the matrix element (0\Q(x) \rj') = m 2 ,f v i / J2Nf. This defines 
the constant and the normalization of the fields in our effective Lagrangian. 
It is interesting to calculate the topological susceptibility 

X v-+oc V 2N f top V (fmi- - \ml)m 2 tov> + 2m\m\ -m%) K ' 

We observe the well-known fact that the topological charge is completely screened if one 
massless quark is present. Note that this formula has been derived for m u = m^- Indeed, 
for zero light quark masses we have = and \ = 0. If m s = we have that = 2m 2 K 
and also x — 0- 

As has been particularly emphasized by Dowrick and McDoughal |20| , we can look at 
the effective Lagrangian in a different way. Namely we integrate over the pseudoscalar 
singlet and are left with a residual instanton interaction. In the present case, with the 
n' — i] taken into account, we have to integrate over both O an d 08- This results in the 
effective lagrangian 

C = Y^YlQiQj ( cos2 <fiD(m w \zi - Zj\) + sm 2 cf)D(m K \zi - Zj\)j . (A.16) 

ij 

The range of the interaction between the pseudoparticles is of order of l/m n . Ignoring 



the the n' — n mixing would have led to the improper conclusion p2| that interaction 
range of 1/mx- Although it is believed that a Yukawa interaction is not strong enough 



|3fl to induce a plasma phase, the relatively long range interaction justifies the study of 



the scenario recently proposed by Samuel [T7| . 
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Figure Captions 

Fig. 1. The integrated topological charge in the Debije cloud of an anti-instanton in a 
slice of length U in the 4 direction of total length 2L. The results are for N c = 3 and two 
flavors with masses (in units of fwT 1 as given in the label of the figure. 

Fig. 2. The pseudoscalar correlator K P (l A ) (upper figure) and the scalar gluonic correlator 
Ks(h) (lower figure) as a function of the length of the slice, U, in the 4 direction of total 
length 2L. The full line is for randomly positioned instantons and the dashed curves 
represent results of a fit. All result in this figure are for two colors and one or two 
massless flavors. 

Fig. 3. The mass dependence of the pseudoscalar topological correlator for three colors 
and two flavors of equal mass (in units of fm~ 1 given in the label of the figures. For 
further explanation we refer the the caption of Fig. 2. 

Fig. 4. The pseudoscalar correlator K P (l A ) (upper figure) and the scalar gluonic correlator 
Ks{h) (lower figure) as a function of the length / 4 of the subvolume for N c = 3, Nj = 3 
and m u = m d = lOMeV and m s = 150MeV. Further explanation can be found in the 
caption of Fig. 2. 

Fig. 5. The mass dependence of the scalar topological correlator, K s {h), for three colors 
and two flavors of equal mass (in units of fmT 1 ) given in the label of the figures. For 
further explanation we refer the the caption of Fig. 2. 



22 



